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Inleiding 

Dit kollege heeft tot doel een Samenvatting te geven van die plasma- 
fysische verschijnselen die voor astrofisici van belang zijn. Weinig 
aandacht zal worden besteed aan randeffekten (van groot belang in 
fusiereaktoren) en veel aandacht aan (energetisch erg onbelangrijke) 
stralingsversc hijnselen. 

We behandelen achtereenvolgens: 

golven in plasma's zowel ín de koude benadering als in warm 
plasma 

stralingsprocessen, niet kollektief en kollektief 

niet-lineaire interakties, zoals koppeling van golven en de versnelling 


van deeltjes door golven. 


De dispersievergelijking 

Verrewee het grootste deel * van het heelal verkeert in de plasma= 
toestand: het bestaat uit volledig geïoniseerd gas. Aangezien de 
onderlinge krachten tussen de deeltjes Foulomb krachten zijn, variërend 
als r°, is de dracht van deze krachten veel gróter dan in een niet 
geïoniseerd medium, dat beheerst wordt door bijv.Van der Waals krachten. 
De deeltjes "kennen elkaar", en werken veel samen, m.a w. kollektieve 
verschijnselen spelen een grote rol. We gaan eerst na welke golven door een 
homogeen magnetoplasma worden gedragen (8 1). Daarna behandelen we de 
straling die door individuele deeltjes wordt uitgezonden en die als 
kollektieve stralingsverschijnselen kunnen worden opgevat (8 2) en 
tenslotte is van belang de interaktie tussenkollektieve toestanden van 


het plasma onderling en met individuele deeltjes 8 3), 


We gaan uit van de Maxwell vergelijking (Gauss eenheden) 


Vv. ES =ánT 1,1 

v. B = 0 1.2 
a> 

mie tl 1.3 
c c, Ct 

WE = -2 BE 1.4 


3 Kd 
Dit zijn de makroskopische vergelijkingen. Vatten we een plasma op als 
5 N > ied 
een verzameling geladen deeltjes: T= E q, ô(r - r;) en 
EE 
: Ee. 
3 = L d; 9, EE Er ERS met vakuum daartussen, dan wordt D= E (e = 1) en 
ï Á 


ER: { u = 1). Vatten we een plasma op als een medium dat polariseerbaar, 
magnetiseecrbaar, en geleidend is, dan geldt 
D= +4 TP af É 1,5 
=$ 
en dan worden bij afwezigheid van uitwendige ladingen en stromen in de 
vergelijkingen 1.1 - 1,4 T = O0, en 3 = 0, We kunnen dus naar believen 
de eigenschappen van het medium onderbrengen in tT en 3 ‚ dan wel in u en € . 
Deze laatste zijn integraal-operatoren en worden alleen in het geval van een 


isotroop, homogeen koud medium tot scalairen. We zullen veel gebruik maken van 


de fouriergetransformeerden van de verschillende grootheden: 





F(k.u) = 4 Fat exp -il.E- wt) r (2, +) 1.7 
en 
P(E,‚t) =d f ek aw erp + iË.r- ue) Mm Ëo) 1,8 
(2m) 


> ie 
Met F (r, t) reéel wordt 


F Kw)t = F(- k, -w) 1.9 
Fouriertransformatie van 1,1 - 1,4 levert dan op, met € = 1, u = Î en 
de plasmaeigerschsvpen in T en 3 8 
É, E &, W) == 4T á TE, ) 1.10 
Kk. B (K,o) =0 Tent 
kxö &, u) == Le Tk, w) - = Ë (k‚w ) 1.12 
kx (E, w) = 5 È K‚w ) 1.13 
Splitsen van T in T Bae a ER Hi in EAN + Ee en vergelijking 


van 1,10 — 1,13 met de fouriergetransformeerden van de makroskopische 


vgl. geeft dan Do )=Ë + ae Tar en via 1.12 en 1,13 
2 
Kx(k zB) +5 D (Ee, w) bni Jet 1.14 


C Cc 








Is Se = 0 dan is 1.14 


2 ' 
kx Kx) +5 E.E=O 1.15 
c e 
E pd B „AL E 
en samen met 1,13 : ( € D. È rade en 1.16 


1.15 is de gol frergelijking in het plasma. Is verder 


dare &, Gj} = 0 (k, w) „É ú, w) : L17 en 
dan is | 
. > Ì 
a È. Bait rl Ege) 1.18 


Met andere woorden: wanneer we weten hoe het plasma reageert op een eleektrisch 
veld (1,17) dan weten we O en dus € en is 1.15 in principe oplosbaar. 


In de praktijk is 


int Es) = (operator g) ÈC E',t') in de gehele ruimte) 1,19 
of wel 
: t 
hit Er Del dr' fat FE, EN ee EE en) 1.19a 


ruimte eo 


Wanneer f ook nog een fuuktie is van E dan hebben we niet lineaire effekten 
(zie 53). Met 1,19a is in algemene zin niets aan te vangen (probeer maar 
fourier te transformeren!) . Dat wordt anders wanneer 
Î.G. Oel a ES dt BE Erste | (Eet) 1.20 
Int gn 


Lan wordt (laat dat zien!) met 1.7 en 1.8 


En K. w) = ol Kw). Ë &, w) 
met 
ej &, w) = f dr dt f (E‚t) exp = i kr -w t) 1.21 


1.20 eist dat f alleen een funktie mag zijn van E-r' ent - t', hetgeen 
homogeniteit in de ruimte én in de tijd veronderstelt; 1.20 beschrijft de 
ruimtelijke dispersie te het medium (via de r-r' afhankelijkheid) en de 
tijdsdispersie via de t - t' afhankelijkheid. Een groot deel van het kollege 
draait om de bepaling van ok, w). ne enen van het plasma worden nu 
gevonden via de oplossingen van 1.15 (gebruik de sommatie over de indices | 
van de komponenten) 

A,, E‚ = 0, met As; = Se (ek; er 6) + es; 1.22 


ij jj, g 


. 


A 








Se 


Oplossingen bestaan alleen Wie det A; = 0: dit levert een vergelijking 
op voor W= W &), de dispersievergelijking. De eigenvektoren van As; geven 
de polarisatierichtingen van È en samen met 1.13 worden B en de polarisatie- 
richtingen van B bepaald, | 

vit 1,3 en 1.4 via inprodukt met Ë resp. met È vindt men de stelling 
van Poynting 
ár 7 „Cixbe-n È 5 ĳ +È, DE) 1,23 

Deze MEER A ‘levert aan de linkerkant de divergentie 
van de Poyntingvektor TE Ëx ä, de elektromagnetische energieflux, en 
rechts de tijdsverandering van de energiedichtheid W. In vakuum is ve 
B + E s Ia een dispersief medium zit de energie niet alleen in de 
fluktuaties' van het elektromagnetische veld, maar ook in de kollektieve 
bewegingen van de deeltjes in het medium. 


Fouriertransformeren van 1.23 leidt tot (zie Stix) 





B {kw =gr Ër xii) ) = Gr Re Ë x Ë) 1,24, 
de flux in de elektromagnetigche energie, terwijl 

tE mel LG tet 1E, 0. Ì 1,25 

15 a he | 

de energieflux in de kollektieve deelt jesbeweging is. E, is het hermietsa 
= == = 4 K == X ak bl : En 
deel van €.( € = Ent Ee» met E €, en Ee, Ee). Is er geen ruimte 
lijke dispersie, Aen is T= 0, 
Verder is 

W [B q ol? energiedichtheid in magnetisch veld 1.26 

B 16 m * fluktuaties ° 
en =p) ú 
(web 
Ì *. 


de snergiealehtheld in elektrische veldfluktuaties én in de kollektieve 

beweging van de plasmadeeltjes, De absorptie van de golven dordk gegeven 

door 
Le ar En et | he 


Deze formules gelden alleen voor zwak gedempte golven, dus wanneer € n klein 


is, Tenslotte is de groepssnelheid v, nog steeds 











6. 
> > 
dk Ë B * “E+P & 
De afleidingen zijn geen sinecure en we verwijzen hiervoor naar Stix en 
Bekefi. 
1.2 Fen paar plasmaeigenschappen. 
1.2.1 De Debye lengte. 
Een geladen deeltje in een plasma heeft een potentiaal veld 
ä . 6 
de + e 1.31 


zodat gelijk geladen deeltjes worden weggeduwd en tegengesteld geladen 
deeltjes worden aangetrokken. Neem een dichtheid van 29 elektronen per 
ee (en een neutraliserende achtergrond van Do protonen En). Voeg een 


„ 


lading eg toe in de oorsprong, dan geldt (met n(r) de elektronenverdeling) 


VP GE) =4 re nl) - Ame ap - 4 te, ô(r) 1.32 
o * 
In zo'n centrale potentiaal geldt volgens Boltzmann dat 
a) = A exp e Ò(E)/ KT 1,33 
met A zó, dat f n (E) dT = N= totale aantal elektronen = no, en dus 
Vv 
is A= n, (waarom?). 1,33 in 1,32 geeft 


2 5 6) > 
Vor) = 4 re no(exp Spe - 1) - áme, 8E) 1,34 
In le orde benadering (e ò(r) << KT) : 


“An ne 
2 > o . > 
PD- P- -4ve, 8E) 1.35 


met als oplossing (1) = eg med dû 








e 
a(r) _ ag = Di EE (-r/)) ‚ met À =À ke jk, 1.36 
e 2 r “deb 2 
ARN 0 Arne 


de Debye lengte. De invloed van een lading is dus ongeveer A deb diep in het 





plasma merkbaar. Bij het aanbrengen van een lading in een plasma treagt 2 
ee 2 KT eas Ie 
zeer snel afscherming op via ladingsverdeling. Zij Vo “EK > en id kaki” 
e e 
dan is À vol ue w is de plasmafrekwentie, en À dep ÍS de afstand die 


deb — 
een elektron in Één plasmaperiode w aflegt. Buiten de Debye lengte zijn 


elektronen nauwelijks van elkaars aanwezigheid op de hoogte, daarbinnen 


wÈl. De elektronen binnen een Debye straal zullen dus nauw samenwerken 








£ q n _ id OE: 
en in kollektieve processen speelt de grootheid No = 0) A deb = Ag Vo lo, 


3/2 —i/2 
Ag 


effekten vecl belangrijker dan in koele dichte plasma's. 


si een grote rol: d.w.z. in hete ijle plasma's zijn kollektieve 


1;2.2 Botsingen, geleidingsvermogen en de Wet van Ohm. 
Botsingen zijn belangrijk, wanneer ze een aanzienlijke impulsverandering voor 


het elektron betekenen: dus wanneer de afstand a tussen een elektron en ion 


2 
zo is dat } mv Ten . De botsingsfrekwentie vv) is dan 
vv) vr ev ps 1.37 
c en 2 É 
e 


De bot singsfrekwentie neefnt dus snel af met toenemende snelheid v‚, m.a.w. 


een snel (of versneld) elektron zal steeds makkelijker versneld worden, Het 
4 nge” wap e? 





impulsverl iesper sekonde is wv crv Wes an = 7 ; wanneer nu een 
elektries veld E wordt aangelegd, á oefent dit een 

2 
versnellende kracht e E uit; is dus E > e B ‚ dan zal het elektron steeds 


v 
minder friktie ondervinden en een "runaway" elektron worden. De kritische 


veldsterkte voor het wegrennen van een gemiddeld thermisch elektron is het 


Zreicer veld E) 


wie -2 
De ke 1.38 
0 


Het ics duidelijk dat de elektronen inude.staart van cen thermische verdeling 
veel eerder weg rennen dan de gemidáelde thermische elektronen. Beschouwen 

. 
we nu een plasmâà als een mengsel van een elektronenvlocistof en een protonen=- 
vloeistof, waarbij we de thermische beweging van de individuele deeltjgs 
verwaarlozen, maar niet de: everiviete stroming van de twee vloeistoffen 


ten opzichte van elkaar. 


Dan geldt: : 


en 














8. 
dv. ER 6 ne, 3 
am test 2% vlee Vp” De eË — Ser, B + 1.40 
waarin P de impulsoverdracht per sekonde van de ionen op de elektronen 


is door de friktie van de twee vloeistoffen. Uit het bovenstaande volgt, 


dat ' 


> = : 
P= vam iv: =v | 1.4 
ce e i e 
Analoge vergelijkingen gelden voor de ionen en verder is 


ò > ÒT + 
TEV. ov =0, mv =d 
waarbij p de dichtheid, tT de ladingedichtheid en v de gereduceerde 


snelheid van het mengsel is. Verwaarlozen van kwadratische termen in de 


snelheid en met mm, >>, (zie Van Kampen en Felderhof): 


Nn 








at : ne 2 
ze =p te Ër xd 1,42 
m e m cm _m, ‚mc m 
e e ed 2e e 
ne ww i 
EN == d Ed € ] hi 0 | 
met en Ë VJ enosE hek 7 5 wordt dic 
e ec c 
dj _e EE EE 
ET on, Vp * Ar % Etevxb 5 n= mers 1.43 
Dit is de gegeneraliseerde wet van Ohm. Voor pe =O enB=0 reduceert dit 
tot 
2 
oÎ __1 2 1, ® + 
ETS mb w, (E - 5 je Je E - LA | 1.44 


In de stationaire toestand volgt de wet van Ohm: Ì = oË; mhd verge= 


lijkingen zijn dus benaderingen waarbij edes Ri ES Mes In een groot 

deel van de plasmafysika zijn we jütee baart Le verschijnselen 
1 

waar —, 


J 


ales 


>> Va» oftewel waar 0 +, Vit 1.44 zalk na fourier analyse 


Le 














9, 
\ 
w à pn 
„ij > Je Ë - v j 1,45 
zodat 
JW) == 7e EON en met 1.17 en 1,18 
w 2 w 2 w Pv 
e(Cu)=l e= rra so I-_E vib 1.46 
c _w w 
1.2.3 Golven in een koud plasma Bo = 0 
Met 1, 22 en 1.46 en vs 0 vinden we 
2 w 
> > w FJ 
kx (k xE) + (1 — B) E = 0 
bad je 2 w 
a) k//E 
Dan is w = Bn Dit zijn de longitudinale plasmagolven in een koud 
plasma : 
b) KL E, dan is u. en 5 + cr? 1,47 


Dit is de En golf transversaal en alleen met w>w p' 


(bereken uit 1.29 de absorptiecoëfficiënt voor deze straling), 


1.2.4 Golven in een "koud" plasma met B, # 0. 
Voor k // Ë en k 1 B blijft voor longitudinale golven alles hetzelfde, 
Voor KiË en k // B: 


)). E=0 1.48 





2 EE 
Dan is weer («2 - er 1 + Ee 
zodat Hi = 0 E, evenwijdig is aan Een l b- 


Verder is de krachtvergelijking 


e& =Ll3xì 1.49 
na ek en 
je 
> i 2 tr 
en dus v x Bo = EN 3 4 1,51 


Dit in 1,43 levert op na fouriertransformatie en Vp. =O 


2 nn : 
w tu B 
iwi = B E-v Ei ss de Î Me H xb 1,52 
ai ikke bn pew 
waarbij b de eenheidsvektor in de richting van 8 is, en 15 de elektron- 


eB_ 


cyclotronfrekwentie te ee 
e 











10, 


zodat 22 
iw “B Ì 
En en 
Eelt stwtne ete Hah tE 1.53 
Pp &Tpe” w 





Inversie van 8! geeft 9 en na invulling in 1,48 de dispersieverge- 


lijking (doe dit zelf; zie Van Kampen en Felderhof): 














Ts w „ Ze ai ie je he ik f 54 5 
w 2 we er wu 5 w 
p P ‚ Pp P 
Of 
ae =] 
cre … Ze, QA 4 5 
w ck” = 8, (1 rr ( a 1) 1,55 
of 
ES: ww? 
u wt uw -Q U, 
e e i 


(Doe de opgaven in Var Karpen eù Felderhof, pag 110). 


2 2 
Voor w 220 A en Q,, 2 > 1, em, golven 
ú u 


Voor w << We Q en Qs wordt 





cr? d w 
in 1 + 5 » US FE aa” . 1,57 
w ni 


en vinden we de Alfvengolf weer terug. 


Voor golven met q, << we IN wordt 1.56 








2 2 

eek- w w 

get MIE ee 
w e e 
Voor deze golven is en 
dw of dd - 

Va me Bk had 7 1,59 

%p 


zodat Je toeneemt bij toenemende frekwentie. Deze golven zijn de zg. 
whistliers, 

In het algemene geval moet men 1,43 fouriertransformeren, voor Vp,” 0, is 
via de krachtvergelijking altijd een tensor 0 te konstrueren die dan de 
dispersievergelijking oplevert. Hierboven is al gebleken dat n = lS geen 
funktie is van k, d.w.z. er i8 geen rimtalitke diopersie. Dit is ook in te 


. k E 4 EA ed . Là & 
zien wanneer men bedenkt dat in een koud plasma in punt x de deeltjes zich wel 


kunnen herinneren wat het electrisch veld was in het verleden (tijdsdispersie), 


hes 








u, 
maar niet weten wat het veld in de buurt was (ruimtelijke dispersie), 
omdat geen plasmadceltjes deed tüformatië kunnen overbrengen; ze bewegen 
immers niet. Zodra thermische beweging van de deeltjes wél wordt meege= 


é ; > 
nomen treedt ruimtelijke dispersie op, en is n wél een funktie van k, 





Zie voor de klassificatie van de verschillende golven in een koud plasma. 
Stix p.9 e‚v. en p. 27 e.v. 


Beginnend bij lage frekwentie vinden we 





l, snelle en langzame mhd golven, met 


2 bmoe® 


4TPC î 
nets ‚ resp, a? coste = 1 + 1,60 
2 B_2 


Bo 0 








2. Bij 9, verdwijnt de langzame Alfvéngolf in de vorm van de ioncyclotron- 
golf 5 de snelle Alfvengolf verandert van karakter, wordt de whistler- 
golf en verdwijnt bij Q. cos 6 ij 

32 In de buurt var Ù beginnen dan de gewone en de buitengewone e.m. 


golf en de longitudinale plasmagolf. 


1.2,5 Golven in een warm plasma. 

In een warm plasma dienen we de stroom in het plasma H te bepalen als 
funktie van cen aangelegd veld Ë, rekening houdend met het feit dat deze 
deeltjes bij het begin niet stil sittenmaar al cen thermische beweging hebben 
(Stix p.l31 e.v.) Dus in plaats van uit te gaan van de krachtvergelijking 
1,49 en de (gegeneral iaeerde vet van Ohm, 1.43, dienen we na te gaan hoe de 
verdelingsfunktie fa (r, v) an de elektronen en f, CF ,v) van de protonen 
verandert. We gaan er voorlopig vanuit dat dit kleine afwijkingen zijn van 


de (genormaliseerde) Maxwell-Boltzmann verdeling: 


a d 3/2 av? 
fv) = GP. exp - Fr | 1.61 


De invloed van elektrische en magnetische velden op f wordt gegeven door 


de Boltzmann vergelijking 








df, df, dq:  & HE, dE, 
a + V Er me ZEE, Ea) 1.62 


Ôr ig & Jv dt C 








12. 


° af, 
waarbij ( me de invloed op fj t.g.v. botsingen aangeeft en i de soort 


deeltjes (elektronen of protonen). 


de ladingsdichtheid is dan 


Te F ayn S Pe d°v 1.63 
en de stroomdichtheid ’ 
Î= Fasn; ff; va 1.64 


Nemen we B, = 0 en O-de orde verdelingsfunkties die onafhankelijk zijn van 
r en Î en verraaflozen we botsingen’ dan geldt voor de le orde funkties ‘f' en 


E' 





p (6) ô Ea 
EN 2 
ee ME ed 1.65 
êr a av 


de gelineeriseerde Vlasov vergelijking zonder magneetveld. 


Fourier transformatie levert op 











ne ij 
-iw PE De 1 ER LEW) -— Sk) - 1,66 
23 m 3% 
Ì Vv 
zodat e Pt a). | 
£'(k‚w,v) = i L es 1:67 
3 m 
> hed 
wk, v 
en á 
et en È &, w) 
Jew) == fa% Es v 1.68 
w = k,v 
zodat met 1 15 en 1,18 
be B, 
ne 
be ese Jee, a en 1.69 
(w -k , v) . 


Dit lijkt heel aardig, maar het-is zinloos zolang we niet weten hoe we 
over de singulariteit bij. w n Kk . v moeten integreren, Hierover is uit- 
gebreide literatuur, zij bijv. Van Kamben en Felderhof en Stix. Zij voor 
het gemak Kk in de z=richting, dan zien we dat voor longitudinale golven | 


de eis geldt 








15. 


òf 
vz 
1+ Pr dv dv dv 5 = 0 1.70 
w XxX y w kv . 


Zonder ons om alle problemen van de komplexe integratie te bekemeren, 
kunnen we nu 1.70 voor een Maxwel1-Boltzmann verdeling voor fo bepalen, 
door EE te ontwikkelen rond v = 0, en zo de Cauchy hoofdwaarde van 
de integraal te berekenen. Verder wordt de pool inw = kv meegenomen in 


het residu: we vinden dan 
_ wp? 3k? ar a MER a 
1-14 ten] an rol exl) 0 171 


En ek 1 
waarin 4% 5 pil Zo 


2vo 


Is ruw zw, + “woz dan volgt voor het reële deel van 1.71 


De 2 2 (KT 
Bmw FIE Gd 1.72 


en voor het iraginaire deel. 
ie) Ln wm 
bend l  Gp  eg 1,75 


'M.a.w. de longitudinale plasmagolven vertonen ruimtelijke dispersie en 
hebben een groepssnelheid 


3k v 
m „ êw 0 


ge dk Be or 1.74 





Zelfs in een "botsingsloos" plasma wordt een plasmgolf gederpt volgens 
1.73. Bezie een plasmagolf met kleine amplitude en fasesnelheid w/k: in 

het stelsel meebewegend met de golf zien we de plasmadeeltjes door de golf 
lopen. Het merendeel heeft een zodanig andere snelheid dan. w/k, dat bij het 
binnerkomen wet energie ren: zen à die bij het verlaten weer vordt 
verloren. Alleen deeltjes met vn E worden in de potentiaalput van de | 


golf gevansen. Aangezien in een termische verdeling zE ‚v altijd < O0, zijn 











uu, 


er iets meer gevaneen deeltjes die langzamer gaan dan E ‚ dan die 

vlugger gaan,Dus wordt door de golf meer arbeid verricht dan ontvangen en 
verd de wolf gedempt. Is echter ergens ze > 0 (een bundel bijv.) dan neemt 
de colf energie toe eri hebben we een bundel plasmainstabiliteit. Deze Landau- 
demping heeft'ook een analogon in het geval van een magneetveld: cyclotron- 
demping (zie Stix hoofdstuk 8). 

Op dezelfde (maar moeizame) manier kan men ook de Vlasov vergelijking mét 
magneetveld hanteren (1.63 met es) = 0}, 

In plaats van 1.68 waar o nog een scalair is, krijgen we een tensor. Uit 

de gelineariseerde Vlasov vergelijking wordt f' gedestilleerd, die verme= 
nigvuldied nk - ve, en „BsIntegreerd over v de stroomdichtheid Ì oplevert, 


en dus 65 mt € jd 5, blijkt dat 


Ma BS en 
Jade Bate A 1.75 
wa i wie 


5 


on M is zelfs voor een thermische verdeling nog zeer ingewikkeld: zie 

Stix p.186 >.v., Si tenko p.9, Belefi p.229, Buiten de genoemde moden 

treden longitudinale moden op bij harmoniesen: van de ion en elektron 
cyelotronfrekwenties: de Bernstein golven. 

Tenslotte treedt in een warm plasma nog de geluidsgolf op: bij zeer lage 
frekwenties is dit een gewone geluidsgolf of ion akoestische golf die bij 

B, # 0 beïnvloed wordt door het magneetveld en een langzame magnetoakoestische 
mode is, bij A, verdwijnt en daarboven vervangen wordt door een modif'ikatie 
van de alfvengolf en can wel elektron akoestische golf genoemd wordt. Hij 


verdwijnt bij Wi . 


1.2.6 Samenvattend hebben we:(zie Kaplan en Tsytovich p.15 e.v.) met dalende w: 





(aannemende dat v SA AS w): 
1) 2 clertromgnetische golven; in principe transversaal, en voor w >> Wp 


zijn dit de eewone lichtgolven met Ve“ e, en v‚ = Cc. 
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w 2 





ext ,e= È,Ë 


WEE) = (— | ) 1.76 
(17 Oef )eos 8 i 
voor quasiloneitudinale propagatie; d.w.z, voor 
w, 
a? sint 0 << wf ( 1- LE. ) cos 6 Lefty 


w 
Dit is haast overal het geval, behalve wanneer k vrijwel 1 B is. De golven 


zijn in de QL tenedering circulair gepolariscerd, en de demping is (met 


ze 3 
No =De ) 
£ U 2 w 2 
vlgs =i WB 1.78 
u) C w2 


2) In de buurt van w worden deze golven anders van karakter en geldt 


<< 
voor f n Wo 


wee w + c? Kk? | 1.79 
E ded ek de 
de fasesnelheid v‚, is dan >> c, en vs en << C, 
kij ne: wr Vp 


3) longitudinale plasmagolven of Langmuirgolven met w v W 


2 2 „2 2 er OAN} ' 
oi so, + 5 kfvo + Asin? 6 1.80 


hd A x ie 2 
Je fes Vn 3 Vo k/w, 1.81 
De botsingsderming y = = v PN ade WIN en de Landau deming 
_ Fre N B =2 
yv=-/n/8 wk Apen? exp — 3(3 + (k Apen) ) 1.82 
zodat voor k «< Ap | de demping snel heel erg klein wordt, terwijl voor 
k > An de demping al gauw groter ís dan an zodat de golf eigenlijk niet 


meer bestaat, 


U) elektron Eernstein golven, met wv ri V,, bestaan alleen voor k 
vrijwel 1 Be De fasesnelheid is ongeveer v‚,, en de groepsshelheid is 
negatief. Voor dispersie verge icingen en demping zie Stix p.223 ov. 
5) de elektromagnetische cyclotrongolf bij As (met Kk // B) die bij 
lagere frekwentie tot whistler wordt en beneden Q, tot snelle Alfvengolf, 


magnetoakoestische golf, 








‚# 
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en Me 
HE ed we A. cos 6 1,83 


Deze golf verdwijnt «dusbij 9, cos Opm later voor w > w‚, weer als e.m, 
golf te verschijrin. 
De whistler golf mot a, cos 9 <w < DN cos 6 is al besproken in 


vergelijking 1.58 e.v. 


IN cos @ 
we RAT TE zl ek, zodat 
Pp 
2 A, cos 8 e?k 
VEV 1.83 
5 f 0 


Uit het beperkte frekwentiegebied volgt verder dat 
1.85 


js ee EA! 
Va < v £ < 20 V A 
De demping gebeurt door botsingen én door Landauderping terwijl in de 


buurt van Q, eyrelotrondemping een grote rol speelt. 


6) da ionakoestische golven, begrmend bij wy; En met w < Oi 
Deze golven worden zeer sterk gedenpt , behalve wanneer ik Nd T,- De reden 
is dat de voortplantinessnelheid ongeveer Ve (er /m,)? is, is Tv Ts 

dan ís jon-Landauderping era sterk. Voor Te >> T, zijn er weinig protonen 


die snel genoeg zijn om de golf te kunnen Landau dempen. Voor Te >> T, is 


v. k 3 vj°k? 5 
wv EE CS CEA + 1. 
(C1+ RT ° pi 


Bij a, treedt slechts ioncyclotron absorptie op, maar voor w < W wordt 
deze golf tot langzame magneto akoestische golf met * 

wv v, kl eos 6 | 1.87 
welke ook sterk gedempt wordt door ion-Landau demping. Nog lager 


nemen botsingen het over en is 
a dee } K: 


7) de jon Rernsteingolven bij n Q, gedragen zich net zo als de elektron 
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Bernstein golven en zijn ook voornamelijk 4 Bj: 


8) de snelle en lengzame alfverngolven. 
De snelle alfvenoolf of magnetoakoestische golf is een voortzetting van de 


whistler mode en daarvoor is 


2 _ BE 
Ee # 04 v 2 sine 1.89 
jef Up 


De demping is bij hoge frekwenties voornamelijk ion-Landau derpine, bij 
lage frekwenties botsingsdemping. 
De langzame alfvengolf bestaat alleen bij w < QD, waar hij door ion- 
eyelotrondermine zwear gedempt wordt. De dispersievergelijking is 

wav Vk eos 9 
en demping is erg gering. . Aangezien de alfvengolf een typische torsiegolf 
is, gedragen door de magnetische spanning en niet door de druk is de Landau 
en botsingsdemping ook erg gering. 
Hiermee zijn de voornaamste golven in een warm plasma besproken, een 


dispersiediagram is gegeven in figuur 1 (sie Kaplen en Tsytovich p.24 en 25). 


2.1 Thermische emissie 


In thermisch everwicht wordt het aantal bosonen met met molfvektor 


k gegeven door: 





NEP) a°k at = (exp BR - 11 a*k aft 2.1 
en 
NE) 0E = Af NR, F) a°® D3 
(2m)® 


Voor hw << KT is 2.1 
N( RE) = HT/ wo En | 2,3 
zodat de energiedichtheid per fasevolume eenheid rer eenheid van volume 


in elke rode wordt 


WELL) Ea = B ak o°r 2.b 
| . (2m)® 
Hieruit blijkt weer dat een plasmamode ís op te vatten als een vrijheids=- 


graad van het plasma waar energie in kan worden opgeslagen, net als in een 
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plasradeeltje. Alleen bij plasmagolven gaat die energie in de elektrische 
en magnetische tuktuaties en kollektieve deeltjesbewegingen zitten, terwijl 
voor een individueel deeltje alle cnergie in de beweging zit: Aangezien 

k< krap voor vrijwel.alle golven, is de energiedichtheid per mode per 
volume-eenheid W < KT Asp en Wein VT, zodat W/W, <N*, mot N het 
aantal deeltjes in de Debye bol. (ca voor verschillende moden en verschil- 
lende astrofysische situaties eens na wat de maximale k-vector is, en vind 


NI). In 2,4 kunnen we d°k transformeren naar dw en de ruimtehoek dî zodat 


2 
Ww,P) date EW fn? GE | an aw a? 2,5 
Bn ec Ur 6) 
met n= ek wi 
Voor een isotroop nlásma is 
2 
Howe r- EE pen 2,6 
2m? e? D 


Deze resultaten zijn ook af te leiden uit de algerene theoric voor fluktuaties 
in een medium met dispersie (zie Sitenko ‚, Bekefi). Aangezien de kwantum- 
mechanische aanpak hier nauwelijks gefundeerd wordt ‚ is de klassieke 
fluktuatiebenadering fysisch gezien juister. De resultaten zijn hetzelfde 


en bijvoorbeeld voor Langpuir golven geldt (per a° 4%) 
: Vv 
Ê 





… BE fg 
A 
KE 1 dun? 
Wasp TZ n° don ê.T 
W, „KE ME 
B B Ve 


waarbij n door de dispersievergelijking wordt gegeven en vp de fasesnelheid 
en Vg de groepssnelheid is. In het isctrcpe geval'is verder 
dk 


WW,P) do = WE) Urk? Ge Wz WE) dw | 2.8 
5 e 
Hel 
en W(k) = ne N (k) ‚ zodat een effektieve temperatuur gegeven wordt 
2m 2 
2m 


Dit tomperatuurbeerip wordt erg vaak gehanteerd om de energiedichtheid 








en 
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in de turbulentie van een bepaalde mode aan te duiden. Bedenken we verder 


met 1.27 dat in het isotrope geval 


ni Be 1 (Ù 2 “we 
an Tr | EK ol ze E49 
dan is voor longitudinale golven (en alleen positieve frekwenties 
beschouwend) 8 
2 
PA … TM 1 Awe ; ” 2 
N, (k) = gn ( An ) [É (%) | 2,11a 


en voor transversale golven, met Wo = n? We vind je 


2 22 A 
N= ( nr ) IË Go)? 2,11b 


2n w° 


waarbij de faktor 2 komt vanwege de 2 polarteatierichtingen van de 
transversale. golven (zie Tsytovich p.80 e.v.). Formule 2,11 geeft het 
verband tussen de É(k) dte we kunnen bepalen uit de golfvergelijkingen of 
op andere klassieke nanier, en de N( k) die hieronder een erote rol zal 
spelen in de vorm van de bosonendichtheid ín een bepaalde mode. We kunnen 
die N(K) het beste beschouwen als een boekhoudkundige methode om bij te 
houden hoeveel boscnen geabsorbeerd resp. uitgezonden worden. In alle 


uitdrukkingen die fysische betekenis hebben valt h steeds weg. 


Spontane en geïnduceerde emissie, absorptie. 
In cen systeem met 2 energienivo's, 1 en 2, met F‚> E) en een stralings- 
veld met fotonen ws; waarvoor E, - E5 = hw geldt volgens Einstein: 


ôn pe; | 
TE == us (1 + ND) + uw NE)np = zE 2.12 


waarbij n, en nj de bezettingen van resp nivo 1 en 2 zijn, N(%) de dichtheid 
geeft van Fotonen die passen bij de overgang 1 »2 en Us D is de overgangs 
waarschijnlijkheid van nivo 1 raar 2. De term U 24 geeft de spontane 
emissie, de term -u‚N do ny de geïnduceerde emissie en de term u, NE)n, 
de absorptie. Voor het stralingsveld N(&) geldt 


ong = ug (44 NÀ) ny-uroN)nzzu (ny (nj -n)NE)) 2,15 


s N << 1, dan hebben we alleen spontane emissie, is N> 1, dan iser = * 


absorptie voor nj > n,‚ en geïnduceerde emissie voor no < n‚, en is 


. 








Ü 


> 
EE eg (y= Ne WN 2.1 





Het enige wat nu nog overblij ft voor het berekenen van emissie-koëfficiënten 
is Us os de overganeswearschijnlijkheid, Verder leggen we er de nadruk op 
dat niet alleen situaties zich voordoen waarbij een deëltje met impuls 


’ 


en energie p overgaat naar toestand p’ onder emissie van een boson ín mode og: 


p' 
B À 


fo} 
In het algemene geval van een deeltjesverdeling £(5) wordt de verandering 


in N . (%) gegeven door 








AN (&) End é APX 
A ef Ee he uy B LED) +, B rk, SEP 2,1la 


waarbij helemaal analoog: aan 2,14, de eerste term de spontane emissie 

regelt, en de tweede term de absorptie en geïnduceerde emissie. We hebben 
De > En 

verder aangenomen, det hk << p, 

Uit 2.1la blijkt det in een gewone Maxwell verdeling (en zelfs in elke 

B A ME heen 

isotrope verdeling) waar ed < 0, geen gestimuleerde emissie optreedt. Een 


òp : 
eis voor gestimuleerde emissie, en dus aanwijzing van het mogelijk bestaan 


van een instabiliteit, is dat ed > O0, hetgeen bijvoorbeeld voorkomt in een 

loss=cone verdelingsfunktie en bij een snelle bundel door een plasma. | 
Aangezien er ook energietoestanden zijn geassocieerd met de diverse 

plasmamoden, kan ook een boson in mode og overgaan in een boson in mode g' 


onder emissie van cen boson în mode o'': 


In het eerste geval is er sprake van spontane emissie door cen deeltje, in 
het tweede sevel door een golfmode. De bijbehorende overgangswaarschijnlijk= 


heden u (p: p', og) en ulo; o', 5") geven meteen ook het omgekeerde proces, 





absorptie en go! fkonpeling 





p' ed Pp en g' > > 0 


TJ ' is 
o 7 Ci o' es Ca & 


maar ook geïnducverde crissie: 





D mk D! en T> — g'. en fol ik => fe Ald 4 ? 
nn | 
pe Ed DS dg af “ae og" g' L > o' 
D, D, Dz 


Het verschil tussen proces A en B is, dat we in A te meken hebben met 

een 2 nivo systeem met een fotonveld, en in B cen proces met 3 fotonenvelden, 
waarbij allerlei vakjes in de k-ruimten kunnen worden bezet. Met andere 
woorden: in B hebben we te maken met 3 meervoudige nivostelsels; in feite, 

zijn in de Keruimte de nivo's kontinu verdeeld. In proces A geldt dat 
| alleen voor het fetonenveld. Zij de overganeswaerschijnlijkheid voor proces 

B: un, duw.s, de kens dat spontane emissie plaats vindt vanuit a°% van mode 0, 


L 


ha à : . « En : . ° 
naar d°k' van mode o' onder emissie in d?k" jn mode oc" (of o'‚k' geruild 





met c''‚k") is | 
> > a > 
Un N/k) a’ ast dk 2,15 


D.w.z. naar analoeie met proces A, en meenemen van de spontane én geïndu- 


ceerde emissie is de kans voor proces B + D, + Dz: 
ds 0) (NE) + 1 Non kt, +1) 4% at af” 2,16 8 
en de kans voor proces C5: en : 
u Noe Er) ME) (NC) + 1) 2,17 


zodat het verschil tussen 2,16 en 2,17 geeft: 


AN (%) , 
8 = Kolk Ne Ml 
mn € tunt Nier Nort DEN Non (N_+1)} 


= Ug ÍN Non = NNT N Non = NJ} 2,8 | 


De eerste term geeft toename van bezetting van o door absorptie, de tweede 
term de geïnduceerde omissie t.g.v. aanwezigheid van Ns in mode c', de 


derde idem v.w.b. a, en de laatste term is de spontane emissie. 








aì, 


Analoog vindt men: 


AN , BN ‚ 
es ent zul torn NN ANN + NO} _ 219 
De a arpumnten È Den Er je Ià zijn hierbij weggelaten, Duidelijk is dat zowel aan 


enersiebehoud als aar impulsbehoud, voldaan roet zijn „d.w.z. e= krak”, en 
ô > . n pn 
wrw! tw, Aanpozien vele kombdinaties kr, k!!, w' en wi' kunnen leiden tot 
Ed . . > Ld 
Kk, w, moet voor de totale verandering in N_() » 2.18 geïntegreerd worden 
>, > pe > > > > ë 
over k' en Kk", or 2,19 resp. over É en È" ‚ en k en k!, Dit levert dan de 


Uinetische vergelijkine" voor de plasmonen op (zie ook Tsytovich ‘hfdst .IV) 





bijvoorbeeld 
LEAK CE Er Pe BON Den CN (Rr 
me Der Dn Uli, PON, CEP)N „l) NN, Cr) 
=NGCON u") — NB) 2,20 
Dn 


+ idem met Ét en Ìr verwisseld, want dat levert ook een bijdrage 

aan dit proces. | 
Behalve de processer: A(gpontane emissie door sen golf) en de daar bijbehorende 
geïnduceorde processen kunnen willekeurige meer ingewikkelde kombinaties 
bedacht worden, Eén belenerijk proces is daarbij verstrooiing 

Bas Ù 

 , st a 
waarbij verondersteld wordt dat de impuls en enersieverandering van het 
deeltje klein is, Dit proces is vergelijkbaar met Do en Dz wanneer we daar 
oa" resp. o' vervangen door p. Zij de overgangswaarschij nlijkheid voor E: 
Up (Pp, o, o') den is de verandering van bosonen Ne, » rekening houdend met 


geïnduceerde emissie (F1) en Absorptie (F2) 


DP > D D= 


> D 
| No en B 
/ 
ot 
F1 Pe 
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a a°%' 


EN 
her Ef wab kN) LB), A) #1) En 


a a°%r 
6 





tr Bk (kt) PP+ ak A EN, ()+1) 2,21 
of, mot p >>hk, be, 
> 
> a 
en zn _& ) fur nk). ze Nr (k') Ee 
geïnduceerde verstrooiing 
+ lun Salt) £ B) EEE 
(2m) 
absorntie van o' ‚of wel emissie van lo} 
= N09) f u, f ®) eet B 
(2m) 
wetens verstrooiing 
7 ' 
En op analoge wijze wordt de vergelijking voor — ge afgeleid, 
Rij de processen A,C, 1…Dj sE F‚ en FP, verandert niet alleen Nd), 
het stralingsveld, maar ook de verdelingsfunktie van de deeltjes £(5) 5 
analoog aan de verendering van de bezetting der nivo's in het 2 nivo 
systeem: vgl. 2.12, 
Tengevolge vern spontane omissie (proces A) verandert £D) 
ND se TE nt, AD 5 
af (an dp 


waarin verondersteld wordt dat hk << p‚ en dat uj maar een zwak verander 
lijke funktie is van D. Val 2,33 beschrijft dan de spontane emissie 
vanuit D naar D= ti, en het terechtkomen in p t.z.v. spontane emissie 

Per al 


uit p + hk, 


Nemen we geïnducoerde processen mee, dan geldt voor de verandering van £5): 


TE 5 dE u(p+ HEK) N KCE +4) — £(D)) 
+ uB,B) NEF (P-4Ì) - £G) 2,24 


Hetgeen na ontwidmeline van u en f in D oplevert, onder weglating van 

















2b, 


E à 
3e machtstermen in k. 5 


dp 
37: 
EP er nb dk Lu BD ak 2) 3 2.25 
JT (2m)? 5 Ap dp 1 _ 


Dit is cen quesilineaire benadering, en het ding heeft de vorm van een 
diffusievergelijking, m.a.w. t.g.v. het stralingsveld treedt er diffusie 
op van de deeltjes in de fase ruimte. 
Op helemaaì analoge wijze vindt men voor het verstrooiingsproces E‚ dat t.g.v. 


spontane emissie: 





det dE get BED 
Ber teek ur B RON, B Altro) ea 2,26 


en t.g.v. de geïnducverde processen: 


AED) _ A ar BNA 
TE EEE NB Net) 


Card 
er as à > > > ap) 
Blk) Cup(p,k,k!) D(K'-E). fj) 2,27 
dp p | 
Wanneer we resp. u en kennen, dan kunnen we de emissie in alle 
pe BE ’ 
e N 


betrokken moden uitrekenen via de vgl. die zE regelen, en de versnelling, 
vertraging of verandering van de verdelingsfunktie f(p) vie de diffusie- 
vergelijkingen 2.25, 2.26 en 2,27. Uit 2.25 is bijvoorbeeld duidelijk, dat - 
=E ‚ hetgeen betekent dat de diffusie 


een afvlakking van de verdelingsfunktie veroorzaakt, 


in het 1-dimensionsle geval za uk 


We moeten nu de verschillende overgangswaarschijntijkheden voor spontane 
emissie door een deeltje, door een golf en bij verstrooiing bepalen. Hiertoe 


moeten we terug raar de klassieke elektrodynamika. 


5. __Stralingsproeesson. 
Voor het berokenen van de spontane emissie van een deeltje maken n 


we gebruik van de vergelijking van Poynting: 


â + 3 
Vi EZ Oei - ne g ' 
5e Vootaal *T VPeotaar JE. aen 
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waarbij Wen P op cen betreffende golfmode slaan, È het elektrisch 

veld van die mode is, en kf het bewegende deeltje voorstelt. Deze laatste 
term geeft de arbeid aan verricht door het deeltje op het stralingsveld 
t.g.v. dat straliresveld. De geleverde emissie is dus 


EE 


=-f ar he 
We vinden de emissie in elke mode o „door het elektrisch veld te ontleden 
in de verschillende mogelijke ‘polarisatierichti. ngen en alleen de komponent 
Ê° behorende bij a te gebruiken. Fouriertransformeren en alleen positieve 


frekwenties nemend geeft 





beed dr sE ad Fo) Fo) + F0). E01 3,5 
Too T (An)* ® 
We lossen nu É in mode g op uit de golfvergelijkine: 
2 h : 
Cs _ -2 8 End, ki 3 E 
Ai E, = (ks! E 8) ge Sa : DE; = Jz (ew) 3. 


Cc is 
waarbij de sometie nie e is gehanteerd ae 3 de fourierontwikkeling 


van de stroom ven het deeltje is. We moeten nu As inverteren, hetgeen 
grote moeilijkhedan oplevert (zie Melrose, Astrophys.& Space Sei.). Aange= 
zien we voornamelijk: in de emissie in een bepaalde mode o zijn geïnteresseerd, 


. . > n se 
veronderstellen we de: polarisatievektoren ©, met; lengte 1 te kennen. Dan is | 








met Ë = ES, E, = Hej. Vermenigvuldiging van 3.4 met e,” geeft 
r… X Te ane lm iw « 
We, A, zer = En (e, jj) 3,5 
en 
2e DE) =2 Rei. eE. 
=lriw (e,j.) 
ze Re((d.85). RE Ek 3,6 
of, met 
„2 2 1 
es Asser KH (ode, ens + EEDE EDE - 12e LT Ket 
+ a at àd ij k? e? 
is 
3 fi Ct HAI 
P = ed dk Re ak 8rriw | (3 ‚, @ )L ae 
ai (2m)* 0 ek wie 


Aangezi. en we uitsluitend in reële strali nz geïnteresseerd zijn, levert 5. 8 








5.1 
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alleen in de polen van de inteerand iets op: dus 


P-tim 1 fati el.) [28 ww) 5.9 


T > eo T TT 0, 
dwte 
( ) 


8 WW (K) 


OE > Fi 5 . za . 
waarbij w(k) de oplossing van de dispersievergelijking c2k?zw?e’ is. 


Pe Kle 


. Fy n 5 ' à 
Aangezien P(X) = Pr u (B, (zie 2,1lla) is met de integrand van 3.9: 
(2m) S 


„Eb - CD en, LED IP 8 wrat) a 


Kw) mT? 





dw? 0 

GG €) 
PD ww (B) 

5.10 geeft het verband tussen onze kwantunmmechrische boekhoudmethode en de 

klassieke clektrodynamische berekening van de emissie Hoëfficiënt. 


Meestal zijn we geïnteresseerd in’ P(w). Duidelijk is, dat 


> > 
UE) u(p,k)m k2aw dn 
P(w,r) dw d2= PI) KA nn eN ike ‚ 3,11: 
La/ak | (2)5 |aw/ak] 


Cerenkov=stralinz 
ZEEENKOV"S LVA LINS 


Voor een rechtlijrie bewegend deeltje is 


J Et) zat SHE) 3.12 
er. 
3 (E,w) z= 2r a va (w-k.v) 5,13 


Voor spiralerendo deeltjes krijgen we 8&-funkties met w-k,v tm , 


waarbij Q de gyrofrekwentie is, Invullen in 3.10 levert op 


ul D, %) 2 E : 1617? aw [W.85? 8 (w-k.v) | z.1l 


( AwZed 
() 


wzw(k) 

Voor e‚m, transversale polven, is de noemer 2w = 25/w ? + c?k? 

Voor longitudinale plasmagolven is de noemer 2 =:2 4/ wo +3 kv? 

Voor een isotroon medium met brekingsindex n‚. is e%= n? ‚zodat de. noemer 


- Àrs _ ck k ke OE en en 
wordt: 2 n u Tos EM W= L= , zodat he de deltafunktie ze kv cos 9. 











| 








NN 


at 


Hieruit volgt dat er alleen Cerenkovstraling mogelijk is in e.m.straling, 
wanneer n >1. 


Laat zien dat in dat geval ES td ki) 








2r2 2 
Peilaf ud (le Se 3 3,15 
et. oe vin? 
n 


Voor plasmagolven, met € = k/k vindt men 





2 
P= f 2mk? dk sine d En v? cos? 8, =kv cos 8) 3,16 
w? a? 
5e e An v/vo 
waarbij v‚, de thermische elektronendichtheid is. 


0 
Reken zelf voor de aardigheid Cerenkovstraling in andere moden uit. 


Met u(p‚k) is ook de absorptie en de geïnûuceerde Cerenkovstraling 


geregeld. 


Remstraling en verstrooiing 

In principe bepaalt de vorm van j& ‚w ) helemal de spontane emissie. 
Rekenen we ki uit voor een elektron dat in de huurt van een proton komt , 
waardoor de elektronbaan iete verandert, dan levert de geperturbeerde 3 
in 5.10 de remstraling op in de verschillende moden, Het komt dus neer 
op het fourieranelyseren van de baan van het deeltje. Veel gereken, zie 


Bekefí , Melrose e‚a, levert op 





8 at ze ik b, 
D =À pi pa _— 
9 w) = hr 3 te N;N, < 7 > fn bi 5.17 


waarbij < El > de waarde van 5 gemiddeld over de elektrondistributie= 
funktie is, en by ‚b, respektievelijk de maximale en minimale botsings- 
parameters zijn. 

Nemen we de stroom t.g.v. een deeltje dat spiraleert in een magneet- 
veld, dan levert fouriertransformatie en invulling in 3,10 de "spontane" 
syn chrotronstraling of cyclotronstraling op, alweer in elke gewenste mode, 
die bepaald wordt door de keuze van € en de oplossing van de dispersie- 


vergelijking. 
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Gaat men de baan van een deeltje na wanneer het in het elektromagne- | | 


tische veld van een golf (wat voor een golf dan ook) terecht komt, dan. kunnen 





we ook die bean fouriertransformeren, en de spontane emissie in een wille- 
keurige mode berekenen. In dat geval hebben we te maken met verstrooiing 
(proces E). Beweegt het deeltje lengzaam, of staat het vrijwel stil dan is 
er sprake van Thomsonverstrooiing en is 

Bnn Bae We 3.18 

dn w 
tin? 

waarin W(w) de energiedichtheid van de invallende golven is. Aangenomen dat 
de elektronsnelheid erg klein is, is de uitgezonden frekwentie van de golf 
ook weer vrijwel w. Is de snelheid van het elektron erg groot, dan is in 
het meebewegend franc van het elektron de frekwentie weer hetzelfde, maar 
in het laboratorium frame is deze Doppler verschoven. Dan is er sprake van N 
het invers Compton effekt, Is de energie van het invallende foton erg hoog, 
dan kan een gedeelte ven de energie overgaan op het deeltje en hebben we het | 
Compton effekt, 


Voor e‚m, straling en het invers Conpton effekt is 


Ser \3 h 
hulk) = en nen Wo) 3,19 
LN w° 


waarbij W(w') de energie dichtheid van de invallende straling is, met 
frekwentie w', en is de uitgezonden frekwentie. Zie voor de diverse u(k) 
van deze processen Kaplen en Tsytovich, en bijv.Bekefi (Tabel 1). 
Niet-lineaire effekten 

Eigenlijk is alleen de Cerenkov straling met recht als een lineair 
proces op te vatten. Maer ook de synchrotronstraling, en endere remstralings- 
processen worden noe: vaak als lineair proces gezien: we "weten" daarbij. inmers 
de exakte vorm van de stroom en rekenen dan de emissie uit. Beschouwen we 
het spiraliseren van de deeltjes in het magnetische veld, dan wel de baans- 
veranderingen t.g.v. de Coulombvelden van de ionen als storingen op de | 


oorspronkelijke baan, dan zijn ook deze processen als niet-lineair op te 
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vatten. Dit geldt noe duidelijker voor het verstrooiingsproces, waar de 
verandering in de stroom t‚e.v. het aanwezige foton wordt bekeken (zie 
Tsytovich p.201 e.v.). Bij de straling gaat het om het produkt d. ee. 

Bij de remstraling is Ki = 3 + EL ‚ met EL de ongestoorde beweging laesies 


+, 
J 


straling) en J' de perturbatie t.g.v. de ionen. EL v E veld van de ionen zo= 


dat de geëmitteerde straling Pv < E?> van de ionen, d.w.z.v energiedicht=- 
heid van het gemiddelde elektrische veld. Bij synct hrotronstralinge is het 


vx 
Ce 


door het elektron gevoelde E-veld v 





v B, zodat j' v B, en Pv Bî, 
dus evenredig met de energiedichtheid van het magneetveld. Ook bij ver= 
strociing is j' v E-veld van het foton, zodat Pv E2 energiedichtheid van 
de invaltende straline,. 

Bezien we nu golf-golf koppeling, dan hoort bij een golf g, een elektrisch 
veld E‚ en dus/een stroomdichtheid j iK Tengevolge van de ook aanwezige golf 
0, wordt Ìá verstoord: js > j 41+ ä rn 5 Ja' 4 BED, zodat de uitgezonden 


straling P Ef EZ , en dus evenredig is met de onergiedichtheden van 


beide golven. Om dt sfhankelijkheid van IEN van Es en E) te vinden moeten we 
dus het niet-lineaire karakter van het geleidingsvermogen nagaan, m.a.w. 
her si 


J pÒ = 5 5) E, + Ë, ) = 60E4 En OE > ta, Ë? +2, E, + É‚? ) + aes 


waarin 5 het lineaire deel van 5 is en 5 het kwadratische, 


1 


Dit niet-lincaire deel ven het geleidingsvermogen zullen we nu behandelen 


in de benadering van een koud twee komponenten plasma. Dan is voor de 


elektronen 
dv > e 1 > 
RRCA EE E +5VxB) 3.20 
en 
FE Vantz=0,dzeny 3.21 


3.20 is niet lineair/en we gebruiken hem om ij als funktie van È, te >... OD 


te testen. Fouriersnalyse van 3, &0 levert op 
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wi) HS Fe) KVP EBU) + ES Fk, H(A 3,22 


dee _ = A37 A37. ent Bk 
Hierin is k= (X‚w) en dA = d ka ko do, do, 8 (ww, Wo). Sk Ke, Ek). 


Hetzelfde woor 5,21 
walk) = kf n(r,) # (k) a | 3.23 
Uit de vergelijking ven Maxwell halen we = Blk) = kx É(k), zodat | 


vlk)x Blk) 
if in vk) x Ú x Èx)) = 


n k Sr 8 1 > > 
=op VO) Ea) — os Bles) le) Ko) 5,24 


Dit in 3.22 levert 


Eed 








En î k 
VK) = a HK) + EE f De Vk) Ee) 


Wm, wr 
Lok > > > ie. 
tf A (ko.vlk,)) (v(k) — mo, Po) 3,25 
Met v(k) = 0) + VK) B werke ‚> volet det de nulde orde snelheid (ín 





een koud plasma!) is 0, en v= en É), zodat: met va in 3.25 ingevuld, 
"e 


w 
de kwadratische toer in E oplevert: 


DE AK 





eN ee af en} ES 5 zalk 
valk) = at Bla) Elr))= zi el arn, Éey) Be) 3,26 


zus . se > > . ke . . DE à . 
Hierin is k = Ky + kos en we zien dat v In de richting k is. Bezien we 


even een transversale colf met É, ES op dan blijkt dat toch de beweging 


> 


E 
van de deeltjes in de richting k, is (als we die golf met zichzelf laten 
interageren), zodat de transversale golf longitudinale oscillaties teweeg 


brengt. Dezelfde procedure loslatend op 3.23 peeft 


ank) = ng F.v, 0) 3.27 
zodat, met de bovengenoemde uitdrukking voor vi 
BAB oe 
n= Se Ì. Bl) | "3,28 
Mm 


Voor transversâle solven is dit 0, maar de volgende orde benadering 


geeft dan 


> 


hee ie 3 (ir al e ee . i AN 
un, (= Of an, (ey )v, eo) + np (ey )va (ko) 5-29 




















31. 
en ereloor mot ä z de + ds TN: 
si be 
JW) sé nov4 (5) = Te É (k) 5.50 
en 
Jo) ze f Ank) Tolke) + nm (kv, (ko) = 
Ze ot OR TR Ae es B 
. Ee 2 > 
En nie En È je Age 4 5 ì 
ze jf dM- AST (ke). Dr E(k)- en Ek) .E(kID} 
m w e 2m “wu, W 
se 1 ee 1eè 
ke 
Symmetrisch maken in 1 on 2 geeft dan voor de niet lineaire stroom 
dichtheid i 5 2 
wo > (k_.E(k_)) 
rd ze en î A k 5 î 2 2 
Jk) 7 Bm. kh a Bae, ).Ëac) id Ea) TW 
€ 12 2 
EE) 
1D) A 1 3.32 


De lastste tm termen komen voort uit 3.28 on geven dus alleen een 
bijdrage voor niet transversale golven. Val. 3.28 levert de gevraagde stroom= 
dichtheid op die ontstaat wanneer twee golven tepolijkertijd voorhanden zijn 
(of de interaktie van cen golf met zichzeif) on deze in de stralingsformule 
ingestopt stelt ons in staat om de oolfkoppeline te berekenen, d.w.z. 
de Up { og, g', of), Dit is verder een hondewerk, en voor de resultaten 


nj 


verwijzen we naar Tsytovich, en Kaplan en Tsytovich. 
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Fia. Í. Qualitative behaviour of the dispersion curves for waves in a plasma which is in a weak 
magnetic field (wo >e- Ope, but va > v‚). The curves are labelled by the superscripts in the 
dispersion relations. 


LN 














Ve 
Woe ‚.k 

Le) 
MN UU Ou  Wue Loe ue Wee Wai 
MiVe Va Vs c Ys C Vr Vre Vri 


Fio. 2. Qualitative behaviour of the relative damping rates of different waves in a plasma. The 
curves are labelled by the superscripts in the dispersion relations. 
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Appendix 
Tables of averaged probabilities for plasma processes 


Tue introduction of emission, scattering, and decay process probabilities enables us to 
find easily the coefficients for emission, reabsorption, diffusive particle acceleration, scatter- 
ing of electromagnetic radiation, and so on. Using these probabilities one can also easily 
write down kinetic equations for plasma waves and particles. In the original papers one 
usually finds complete expressions for these probabilities as functions of many variables. 


The use of such expressions in actual astrophysical applications requires additional calcu- 


lations which are not always accessible to, say, an astrophysicist or experimentalist. 

In the tables which follow we give an as far as possible complete summary of averaged 
expressions for the probabilities of those plasma processes which are most important under 
astrophysical conditions. A large part of these expressions has been evaluated anew. 
However, even in those cases where the probability for the processes was known we have 
averaged over those parameters, the dependence on which is unimportant for applications. 
For instance, when evaluating the non-linear scattering of waves in isotropic plasma 

turbulence it is sufficient to know the probability averaged over the angle between the 

particle velocity vector and the plasma wave wavevector. Besides, even for the case of 
anisotropic turbulence we can fully ùse such an averaged probability as there is normally 
po particular value in calculating factors of order unity in astrophysical applications. 


„ 


The averaged expressions for the probabilities are much simpler and physically trans-_ 


_ Jucent and can be used to obtain estimates in the overwhelming majority of cases of analy- 
zing plasma processes in astrophysics. We have described in this book in detail the method 
for using the probabilities for various processes. vj 

In the following tables we give expressions for the averaged probabilities ù (in some cases 
retaining also the non-averaged expressions) in the notation common to the whole of the 
book. For each process we indicate the limits of applicability of the expressions obtained. 


|. Probabilities for the emission of plasma waves 


We give the probabilities for the emission of waves of type o, averaged over the angle 
between the velocity vector of the emitting particle, v, and the wavevector, k, of the emitted 


‚ wave. 
For the case of isotropically distributed particles with a given momentum distribution 


function f, these probabilities enable us to determine at once the emission coefficient per 
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unit wavenumber, J2 and the damping rate y°(k) by a simple integration over the absolute 
magnitudes, p, of the momentum. We have 
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gee ij Elk, fs dp, EN 

5 t d 

MB fa nd (bar lle | AN 
wk) = Í u (k, v) dp (5) P dp KAD à Ik (k‚v).fp dp, (A2) win E 


where T is the temperature and the last equality is obtained for a Maxwell distribution. 

One can obtain exact values of JZ and y°(k) using the averaged probabilities #(k,v) …— 
only for isotropically distributed particles, but for obtaining estimates these quantities can 
be used also for arbitrary distributions. In particular, we find for the beam instability growth 
rate: EE: 
‘ f_ Fo*/vove) | | 

Smile} ee 


where v, and Áv‚ are the average velocity and the spread in particle velocities in the beam 

with density n‚ and mass m,. For relativistic beams with a large energy spread the quantity 

m‚{Av)* must be replaced by the average energy €, ha 
We can similarly afso evaluate other parameters. For instance, we find for the coefficient 

of diffusive Acceleration in a strictly isotropic plasma turbulence 


De = [52 WILOO (A4) 


where WZ is the spectral energy density of the plasma waves. 
In Table I we have given the probabilities for Cherenkov emission of all plasma waves 
(except electromagnetic waves) for the case when the wavelength is much smaller than the 
Larmor radius, that is, k >> w/v — eHc/ev (unmagnetized particles) since only in that case 
is tbe Cherenkov emission probability not small. In what follows we must thus bear in 
mind this lower limit of the wavenumber for the applicability of the expressions given here. 





: TABLE 1. Probabilities for the emission of plasma waves : 
3 e s É : Ù 
Plasma wave mode Averaged probability he 
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1. Longitudinal plasmons (w' = 








me Wi í 1, @plk <= v; 
0, 











Ed 

@‚,) in a weak magnetic field zl = 

(ze > Wz). . 2nv k @y/k > v. 

“2. Gyro-frequency _plasmons EN oke 1, Er Ogelk Sv; 
(wt = wa. |cos |) in a weak magne- = 5 ( v ) 

tic field (pe > or.) bns & W)t Walen 
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en 
Plasma wave mode Averaged probability 


Maximum probability fork = wle, 











vs CU zel2se nv \ oe 
3. Longitudinal plasmons (wo! = gn 1, Wgelk < vs 
ze 09. [cos D|) in a strong magnetic ä = FE Me Ope kv \t 
feld (ao, « We). 4Anv k \ oe) ‚ -Oalk =v. 
Maximum probability for non- k 
relativistic electrons with k S wael?- Ot met (Oe 
For ions and relativistic electrons Uaax “4 ri ( zt) 


the value of äls, is (e/mc®)* times 
larger. 


4. Gyro-frequency plasmons 








































_ 2e me $ of, 1, zelk 8 
Er ee 
field (ze & Og). he e ú e 
kv \? ke\* opt 
5, Whistlers (Ww” = Wz k°X a In (5) +Gin 2-25) . SEE <%; 
x|cos 1 /w?,) in a weak magnetic Die ie ee „ a 
field (oe >> 0). nv k ( ke Ni oiv y wak, 
A ge: get k, 7 We é 
6. Ton-sound waves (@'= zm @ {1 2 ri <v; 
2 1/2 H ze ) { FR MiVs { H 2 % “Vs 
je me daj: at frequencies | #5 vi \0 “Tie A0 if keld), nee 
ate ER Dg 
7. Alfvén waves (* =s vak {cos #1) 
in a magnetic field satisfying the ze mi TAR 1, VA SU; 
condition c > va >> v‚. The condi- Pe oe (2e) valk +20) J 4? B 
tion for the applicability of the 5 pa 
_ equations is oafv < k < vole. 
Maximum probability for particles zr _r mvv’ (22) 
with Vv Ui © mt 8 nc \op/° 
8. Fast magneto-sound waves 
(w* =— v‚k) in a magnetic field satis-" zm (0a\? 1, va <9; 
fyingtheconditionc ze v4 > Vs: SE (Ze) vvv) Le ape k 
The limits of applicability of the 5 ke 
formulae are: wg/v < k < wyjlc. » 
9. Slow magneto-sound waves 
(o®* = v‚kleos df) in a magnetic 1 u pt 
« E e , . 3 Ed 
field satisfying the condition mn FMS 
CVV, The limits of appli- EFA ne (5) WENT oet 
cability of the formulae are: Vv, Wd 
oplev <= k sono. ° Hú : 
es … Mies Te ú 
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Appendix 


Mi. Probabilities for non-linear scattering of plasma waves by epithermal particles 
leading to the generation of electromagnetic waves 


“ When Langmuir and gyro-frequency plasmons, and also ion-sound waves, are scattered — 
by epithermal particles high-frequency radiation is produced with frequencies w >e- wo, Of 
@ >> Oje The particles are assumed to be unmagnetized, v => w/k. We give in Table II 
the probabilities averaged over all angular variables, that is, both over the angle between … 
the velocity vector and the waveveetor of the scattered plasma wave, and over the angle 
bétween the velocity vector and the wavevector of the emitted electromagnetic wave. In the 
case of isotropic distributions of the particle velocities v and of the plasma wave wavevectors 
k, the emission coefficient per unit frequency can be evaluated using the formula 


: Ee 3 fr zi me: 5 dk: dp 
| JE = WE [ dw, ka, ofoWi, oki) : (A5) 
ô ò Ë 


We have used here for the dispersion relation for the electromagnetic waves w == ck; kr 
is the wavenumber of the plasma waves with spectral energy density Wz, 
__ TABLE út. Probabilities for non-linear scattering of plasma waves by epithermal particles 


Conversion process Averaged. probability 








IE. electrons (l-+e == t+e”) 
1. Langmuir wave scattering ina | , a @… (@0,.\? @ wt 1, oki < vas 
weak magnetic field (@' ss w,). | à = do 5 ( pe ) rd Ps (oto) ha al, ie 


In a strong magnetic feld (w' =s ad 











€), Icos DI) the conversion proba- Oe dons (li t+i) 
bility is smaller by a factor two. ’ En TE Ope (ze y Ope 1, w/k, < 03 

6 ni \o/ kwe, O, w/k, — vj. 

electrons (he =t+e 
2. Gyro-frequency plasmon scat- ì Mt ( 8 ) k en 

tering (w* = op, |cos 9) in a weak | zt — Z. Ppe (22e) Ope (eo) 1, ok, < ves 
magnetic field (@,, 2e 07,). In the 5 a le / kweet ki 0, w/k, > ve- 
case of a strong magnetic field |- -_ fons (h+i == t+i/) 








(@? = op) the probability must be úa ië ore ee) on { 1, lk, < os 


Ri s _ Ope 
multiplied by (w‚/os)“. dn = dr Ea Ôreilkneen: 


electrons (ste == t+e’) 


ES TE Ope (5) Mi pe ( 8 ee) 1, ofki < vs 
Ds = oe vi 











3. Scattering of ion-sound waves | __ 10 me \ol me kwe \ Ki) LO, o/k,= ve. 
_@= vka} + kid) HW). Á ions (s+i == t+f’) 
zt TH Ope (5 ) MM; pe { 1, o/k; — vis 
kT 6 a \e/f me kw 
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Similarly we have for the damping rate (due to the inverse conversion of electromagnetic 
waves in plasma waves of type 0): 


íC jo dkad 
Ze | Í Beke Pig eer | | Vetere: 
0 © 


(A6) 





Sinilarty we tan obtain approximate formulae also for anisotropic distributions. 


V. Coefficients for the differential transfer of plasma waves within one mode or ie 
between different modes through non-linear scattering by thermal ions and age di 
fusion processes 





Averaging the scattering probabilities over the angles and assuming all plasma tuttle 
‚modes to be isotropic we can obtain the differential expressions to determine the change 
with time of the spectral density of plasma waves for the case of differential transfers. 


TABLE v. Coeficients for differential transfer of plasma waves 





Formula for the change in the Scattering and conversion 


Scattering or fusion processes 














spectral density coefficient 
all = al = 
. me 
e 
1. Transfer of longitudinal plas- 1 27 nemo 
mons (@! =: 0.) within one mode Di dS om E ware eet 
Atie Ii, Hs = 1. ‚öt ek 1 mir 
fi U#7/TD* Ld . De | 
. t, …_ TT. 
2. Transfer of longitudinal plas- 
mons (@' <= w,) into transverse ai PW? Li 5 
plasmons (or 3 je) (Hi == p'+Y, Dt Wi — oz CWD. air —= at. 
1+4s = p’). 
3. Transfer of transverse into lon- 1 B 
gitudinal plasmons (pis I+4i, 0: Wk DE (55) ‚ Bi: Eea 
pts == |I). ôt k 1’ aP} = oe}, 
282 








Appendix 


OG Ons 


Scattering or fusion processes 


ET OO 


4. Transfer of transverse into 
transverse plasmons (p+i ze p’+iï’, 
„Pts =p). 


S. Transfer of high-frequency 
transverse waves (@ zo 00) (LHE ze 
zette’, ttl = 0). 


_ 6. Transfer of whistlers into 
whistiers (wi == W+Ï). 


7. Transfer of whistlers into ion- 
sound waves (w+i == S+Ï’). 


8. Transfer of ion-sound waves 
“_Înto whistlers (si == w+i’). 


9. Transfer of ion-sound into 
ion-sound waves (s+-i == 8’ +’) and 
of slow ‘magneto-sound waves into 
slow magneto-sound waves (ms +i == 
ze MS + i°), 


10. Transfer of fast magneto-sound 
‚and Alfvén waves both within one 
mode and from one mode to another 

‚(aki zr ati’; miss mi’; 
mii ze a’ +i’). 


11. Transfer between the Alfvén, 
fast and slow magneto-sound wave 
modes 
(mi z- msi’) 

(ati se msi’) 
(mshi z= m+i’) 
(msti ze ati’) 


Formula for the change in the 
spectral density 
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Scattering and conversion 
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When evaluating the transfer between transverse and longitudinal plasma waves we have 
used in Table V the “longitudinal” wavenumber scale (see $ 3). 
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VI. Coefficients for fusion and decay of waves, leading to the generation and absorp- 
tion of electromagnetic radiation near the frequencies w‚, 20 and 2m 


Averaging the fusion and decay probabilities over the angles and assuming all plasma 
turbutence modes to be isotropic we can obtain expressions to determine the change in time 
of the spectral densities of plasma waves in fusion and decay processes of different modes 
(Table VI). 5 


TABLE vr. Coefficients for the generation of radiation through the fusion of plasma waves 





Fusion coefficient 


Eusion and decay processes Formula for the change in | 


spectral density 








1. Eusion of longitudinal ptas- 
mons into an electromagnetic wave 
































with frequency wt ss 2wll+Y — t). OW Wa, \* dk ‚…_ Ív3z Wie 
The phase velocities of the longi-| - ar “PJ \E, » EE 
tudinal waves are smaller than the i 
velocity of light. 
2. Fusion of two gyro-frequency ie ae 
plasmons into an electromagnetic ow = ft í (2) dk pr = 322 oere 
wave with frequency wt S 2wor.lh+ ot ks zc hl 5 none 
h’ hae ĳ. st « 
3, Fusion of a longitudinal plas- P zel ú 2 
mon and an ion-sound wave ONE ge MME Be = ES 
(ls — p). ôt k 2V3 namvret 
4. Decay of a transverse plasmon,{ ____‚‚p pis 3 ke 
involving an ion-sound wave ow: = pre WaW: Ni Br = EE ED pal: 
BES). 7 ot k_ | 8 namdre 
ee ee Sne 
5. Fuston of a longitudinal plas- 8 Wi =p WiWs B = v2 em 
mon and a whistler (l+w — p): Ot Ke oik name” 
6. Decay of a transverse plasmon, | _ OWZ „ WiWi EDE” 
involving a whistler (p — |+ w). : or B" ks pr Ei 8 oberemn e 4 
7. Fusion of a longitudinal plas- we? ra 1 at 
mon and a fast magneto-sound or = ge B 5 pr = — Ss: 
wave (1m —> p). t V2 hemiva-c 
8. Decay of atransverse plasmon owE z où 
involving a magneto-sound wave ern =P" WWE, pen = ri er 


(Pp — 1m). 
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Vil. Probabilities for fusion and decay of high-frequency electromagnetic waves 
with plasma waves 


We give in Table VII the probabilities for the fusion and decay of electromagnetic waves 
with frequencies w ze @pe and w => @yje involving all possible plasma turbulence waves 
with wavevector k,. Only the coefficients are averaged while the delta-functions take into 
account the conservation laws during the fusion and decay processes so that the probability 
Up, w‚n, fOr the process can be expressed as follows in terms of the quantities (wo, k‚) given 
in the table: ’ 


ue won, = (2) ala, ka) A Kk) òke) — 0'(l")—0(ke1)). (A12) 


One must substitute expression (A.12) into equations (5.20) to (5.22) after which one can 
in particular, obtain the subsequent formulae from $5. 
To determine the cross-section for the elastic scattering of electromagnetic waves with 


o => ck, one can use the simple approximate formula p 
! 2x ck 
0 = en Í a(w, ki) Wi, (ad (A.13) 
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Tante vu. Prababilities for the decay of high-frequency electromagnetic waves _ 


Sne 


Fusion and decay processes 





Averaged probability coefficients 


1. Fusion and decay of electromagnetic waves 1 & okt _ 1 oek? (Ope \* 
involving longitudinal plasmons (t == t+I). en ne Tr ( ) . 


In the case of a process in a strong magnetic 
field, where w' = w‚‚lcosd|, there appears a factor 


1 & wilt 


18z mi oz 


2, Fusion and decay involving gyro-frequency 
plasmons in a strong magnetic field (t == t'+h). 





1 & oki 





3. Fusion and decay involving whistlers 


„ 


4, Fusion and decay involving ion-sound waves 





(t == t’+5). 12 mi vet 
S, Fusion and decay involving Alfvén waves 1 & ck 

(t == t+a). ik 48 mi vo? 
6. Fusion and decay involving fast magneto- 1 & wik, 


sound waves (t == t’+m). 
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